DSF-44/2003 
CERN-TH-2004-015 



Bounds for neutrinoless double beta decay in SO (10) inspired 

see-saw models 

F. Buccella and D. Falcone 

Dipartimento di Scienze Fisiche, Universitd di Napoli, 
Via Cintia, Napoli, Italy, and INFN, Sezione di Napoli 

By requiring the lower limit for the lightest right-handed neutrino mass, obtained 
in the baryogenesis from leptogenesis scenario, and a Dirac neutrino mass matrix 
similar to the up-quark mass matrix we predict small values for the mass and 
for the matrix element niee responsible of the neutrinoless double beta decay, rriu^ 
around 5 • 10^^ eV and niee smaller than 10~^ eV, respectively. The allowed range 
for the mass of the heaviest right-handed neutrino is centered around the value of 
the scale of B - L breaking in the SO(IO) gauge theory with Pati-Salam intermediate 
symmetry. 
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The robust experimental evidence P, 01 for neutrino oscillations , with increasing 
precision on the square mass differences and mixing angles 

Ami = 7 ■IQ-^eV'^ (1) 

(tan^,)2 = 0.4 (2) 

Ami = 2-6 • IQ-'^eV^ (3) 

(tan^,)2 = l (4) 

should be a milestone for the extension of the standard model. Indeed the smallness of 
neutrino masses points in the direction of theories like the unified SO(IO) model j4|, where 
the see- saw formula jsf 

rriL = -m^M]^^mD (5) 

arises naturally with the elements of the Dirac mass matrix of the order of the charged 
fermion masses and the elements of the Majorana mass matrix at the scale of the spon- 
taneous symmetry breaking of B — L. Right-handed neutrinos of high mass have been 
advocated in the baryogenesis from leptogenesis scenario ^ and a lower limit for the mass 
of the lightest neutrino around 5 ■ 10^ GeV has been found 

In a recent paper Akhmedov, Frigerio and Smirnov studied the compatibility of 
the seesaw mechanism with what is known on neutrino oscillations and the leptogenesis 
scenario, and reached the conclusion that the two lightest states should be almost de- 
generate. Here we study the consistency of the see-saw formula with the leptogenesis 
framework within the hypothesis that m^ has strong similarities with m^, as expected in 
SO (10) theories. We do not require strict equality (which would follow from the assump- 
tion that the electroweak Higgs transforms as 10 representations), since it does not hold 
for md and me, but assume, as in P, the eigenvalues of mn to be 1 MeV, 400 MeV and 
100 GeV, respectively, with the same hierarchical pattern as the up quarks, and small 
mixing angles for the matrices Vr and Vl which diagonalize m^ in the basis where me is 
diagonal: 

mn = Vr diag(mz)) V^- (6) 

We define 

Mh = {diag{mD))hh (7) 

and mi (i = 1,2,3) the eigenvalues of mi- One has for the mass of the lightest neutrino 
the upper limit 

, , , (MiMaMg)^ ^ „ 

< = 'jfkj^ < (0.23 eVf (8) 
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comparable with the hmit coming from astrophysics P. From Eq.(6) and the inverse 
see- saw formula 

Mr = -rriDml^m^ (9) 

it is easy to derive 

M5, = -V,'}V,^MfM,A% (10) 

where 

A% = {ml'UV:/^V,f (11) 

which shows the intriguing property that the part proportional to MfMg of all the matrix 
elements of Mr is also proportional to the same factor A^^. This part may be diagonalized 
by the transformation induced by any unitary matrix with (Vff , Vgf , VJjf ) in the last line. 

If we take Oatm = 45°, 6^13 = 0° and define c = cos 6^5, s = sin6's, we get the following 
elements of the effective neutrino mass matrix: 

rriee = c^mi + s^m2 (12) 



rUe^ = cs{mi — / \/2 (13) 
rrier = cs{mi — 1112) / \/2 (14) 



= {s^mi + c^m2 + ^3)72 (15) 



^iiT = {s^mi + c^m2 - ni^)/2 (16) 
nirr = {s^^i + (?m2 + m^)/2, (17) 

and a similar expression for m^^ with the exchange rrii l/rrii. From the expression for 

we get 

^33 = [cvi* + -^{V2f + vi*)?:^ + [sv,i* - -^{V2f + vi*)?:^ + l[v^3 - vir^ = 

y I "ii v2 '''■2 ^ ''''3 

^[v^3r]'(— + —+—)+••• (18) 

/ m\ 1712 ^3 

where we have selected the term proportional to {V^*y, since we assume that, as it 
happens in the CKM matrix the largest matrix elements are the diagonal ones. 
From now on we shall take tan^ Og = 0.4 0. 

If \mi\ ^ (tan^ ^5)1^2!; the lightest eigenvalue of Mr would be M^ ~ M|/7mi and 
then the product 

M^M^ ~ '^[^^^^y ^ 2.6 ■ 10^5 GeV^, (19) 
m2m3 
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which imphes < 5.1 • 10^ GeV, an order of magnitude smaller than the lower limit 
required in the leptogenesis scenario. 

To get some cancellation in the r.h.s. of Eq.(18), \mi\ should be at least ~ (2/5)|m2|. 
This would imply for the product of the eigenvalues of Mr 

MiM^M^ < 21(^3M2Mi)^ ^ lO^OGeV^, (20) 

and therefore, in order to be > 5 • 10^ GeV, one should have 

Mg^ < 4 • lO^^GeF (21) 
which requires, according to Eq.(18), 

^2 2 

I— + — + — I <0.8eV^-^ (22) 
mi m2 

and a cancellation between the three contributions, in fact > 4 eV~^. We should also 
have 

M^M^ < 2 ■ lO^^GeV^ (23) 

which implies that the coefficient of (M3M2)^ in the expression of M^M^ be less than 
1.25 cV-2. 

For hermitean matrices, if the lightest eigenvalue is smaller than the other two, the 
product of them is approximately the sum of the minors corresponding to the three 
diagonal matrix elements. Let us consider then 

(24) 

where A^^ = Vj^V^^m^j, which shows that the terms proportional to (M3M2)^ for all the 
minors considered are multiplied by the factor 

' ^ ■[c^n + 4(^2i + ^3i)]'+ 



[{Vi-Vi)f\<1.25eV-'. (25) 



Since we have 



2mim2 

1 



> 25 eV-^ (26) 

777.2777.31 
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to get the coefficient of (M3M2)^ smaller than 1.25 eV~^, we need a cancellation between 
the three terms in Eq.(25). Let us see if it is possible to obey at the same time the limits 
given by Eqs.(22) and (25) with given by Eqs. (13-18) and the form 



^ cos p sin p ^ 



sin p cos p 

V ly 

which looks like the CKM matrix up the order A. Eq.(25) takes the form 
I — ^ — (c^(cosp)^ + s^(sinp)^ — a/2cs cospsinp) + 
1 

(s^(cos pY + —(sin pY + \/2cs cosp sin p) + 



(27) 



mim2 2 
1 (sin p) 
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, < 1.25 (eV)-^ (28) 

171x1712 2 

In |ll| one found the solutions of Eq.(22) with a vanishing r.h.s. consistent with 

ml-ml = Ami = 7 ■ 10"^ eV^ (29) 
ml - (?m\ - s^m\ = 2.6 ■ 10"^ eV^ (30) 



one with the same sign for m2 and ms and 
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!:^i = -(ta„«.y^-^.-3^. (31) 

777.2 m^C^ & + U.lo 

the other with the same sign for m\ and ms, 7773 larger than the other two, and 

777.1 — — ^^2, (32) 
3 

7r7i ~ cos20s^3 = ^"T-S- (33) 

These two solutions are shghtly modified by the small upper bound in the r.h.s. of Eq.(22) 
into 

7771 . ^ ^^c^ ±^.%- eV-^J Ami -2 ±0.02 

— ~-tan0,^ =^ - = 34 

7772 ^ ^ ^2 , 6.12 ^ ^ 



mi 



cos 2^,(1 ± 0.8 ■ ey-V^^a) = -(1 ± 0-04). (35) 
7773 7 

There are also solutions, with 777.1 and 7770 with the same sign and |r77i| > 0.11 eV or > 2.5 

eV (which is excluded by astrophysics p). However only the solution of Eq.(22) with the 
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same sign for m2 and ma satisfies Eq.(28) with a small value of p. By considering only 
the term proportional to {V^Y in Eq.(19), one neglects terms at the order 

(36) 

for the relative coefficients of — and — and 

mi m2 

% (37) 

^33 

for the relative coefficients of — and — . Should one take for the ratios of the matrix 
elements of the corresponding values of the CKM matrix, these uncertainties would 
be 10% and 100% and would effect the first and the second term in the denominator of 
the r.h.s. of Eq.(31), respectively. By taking into account all the uncertainties one has 

- 0.45 < — < -0.25 (38) 

which implies 

4.23 • 10"^ < mi < -2.16 • 10"^ eV (39) 

8.64 • 10-^ < m2 < 9.4 • 10"=^ eV (40) 

5.15 • 10-^ < ms < 5.18 • 10"^ eV . (41) 
We predict a small value for the neutrinoless double beta decay parameter 

Mee = sin^ ^sm2 + cos^ Ogmi — — 17111712/1113 ± 0.8eV~^mim2 (42) 

which implies M^e to be in the range (3.5 — 8) ■ lO^^eV and a slightly larger value for the 
effective electron neutrino mass 

m^^ — sin^ ^s|m2| + cos^ ^s|mi| (43) 

in the range (4 — 5.6) ■ lO^'^eV. For simplicity we have considered real m^'s, but our results 
would not change by considering complex m^'s. 
By taking ^433 = and A^^ = one gets 

V2sc{m2 - mi) 

tanp = — ^ ' 44 

ma + c^m2 + s^mi 

in the range (—0.15, —0.12) about twice smaller than tan^c- 
^Prom the identity 

^33^22 - (^23)' = iii/mim2m3 (45) 
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it follows that, if we ask that no term proportional to M| is present in Mr and no term 
proportional to is contained in the diagonal minors of M/j, there is no term 

proportional to (M3M2)^ in M/j. In conclusion, the largest contribution to Mji would be 
the ones proportional to M3M1 and M| which are multiplied by the factors A^^ and A22, 
respectively: 

M^, = -MsMiKfnMi - M^VJnf ^22 + - (46) 
where we have not written the terms proportional to M2M1 and M^. In general 

DetM^ = -(MgMaMO^DetA^ (47) 

and therefore 

DetA^ = (mimams)"^ (48) 

For real M^, which is not what we want to get the CP violation needed to produce 
baryogenesis, but is a property of the CKM matrix at the first order in sin^c, we have 
the eigenvalues of Mr as —M| and ±(M3Mi)/ a/ A22^3m'2"^i- Since 

sinV(— + — ), (49) 

mi 777,2 

the eigenvalues are in the ranges (5 — 6.8) ■ 10® GeV and ±(4 — 5) ■ 10^'' GeV, respectively. 
The lightest right-handed neutrino is just above the lower limit 5 • 10® GeV from leptoge- 
nesis. We have deduced the consequences of = and A^^ = to display an example 
with the lightest right-handed neutrino mass larger than 5 • 10® GeV. 

In general the mass of the heaviest right-handed neutrino will be in the range 10^" 
- 5 ■ 10^^ GeV with the geometrical mean of the two extrema very near to the scale of 
B — L symmetry breaking in the 5*0(10) unified theory with Pati-Salam 0| intermediate 
symmetry, 3 ■ 10^^ GeV jl4'. 

Our assumption of small mixing values for and on the value of the product of the 
two smallest eigenvalues of ( 4 ■ lO'^MeV'^ ) excludes the high values of proposed 
in the hterature [l^. 

The upper limit given by Eq.(21) implies a small coefficient for the contribution pro- 
portional to M| to the matrix elements of M^, which is welcome to get the lepton 
asymmetry needed for leptogenesis, which would vanish, should one neglect the other 
contributions. 

We want to stress that the small values found for A33 and A^^ are not unnatural 
tunings; just the opposite: since Eq.(5) gives rriL in terms of (and of m^), the fact 
that the matrix elements of ( as well as its eigenvalues ) are not so different, as it 
happens for rriL, makes the choice described here a natural one. 
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